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Abstract. We generally study the density of eigenvalues in unitary ensembles 
of random matrices from the recurrence coefficients with regularly varying 
conditions for the orthogonal polynomials. First we calculate directly the 
moments of the density. Then, by studying some deformation of the moments, 
we get a family of differential equations of first order which the densities satisfy 
(see Theorem 1.2), and give the densities by solving them. Further, we prove 
that the density is invariant after the polynomial perturbation of the weight 
function (see Theorem 1.5). 



1. Introduction and main results 

Statistical behavior of eigenvalues of random matrices was first studied by E.P. 
Wigner in 1950s to get an information about spectra of heavy nuclei, and the famous 
semicircle-law was found then. In 1960s random matrices were intensively developed 
by E.P. Wigner, F.J. Dyson, M.L. Mehta and others for a better understanding of 
statistical behavior of energy levels in nuclear Physics, see [Po] as a collection of 
early papers. Later more and more importance was gained in other areas of Physics 
and Mathematics (see [Me] for a comprehensive introduction, or J. Phys. A 36, 
2003). 

It is of special importance to observe that the eigenvalue distribution of unitary 
invariant ensembles which can be described in [FK] by weight functions from the 
families of classical orthogonal polynomials. Such ensembles can be defined by the 
probability distribution on a space Hj^ of A'^-order Hermitian matrices as 

N 

P{dA) = F{A) dajj W ddiajkd^ajk (1-1) 

where F is an intcgrable nonnegative class function on H^. 
Using Proposition 2 in [HZ] we obtain 

/ F{A)dA = 

•••/ i^(diag(2;i, . . . ,xjv)) ]^ {xj ~ Xk)"^ dxi . . . dxN (1-2) 



z 



N 



l<j<k<N 
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where 



^fe(fe-l)/2 



Especially when 

JV 

F(diag(xi, . . . ,a;jv)) = Z'"" JJ ijj{xj)xi{xj), (1.4) 
one introduces the joint distribution function for the eigenvalues 

-PAr(a;i,a;2,...,a;Ar) = ]^a;(a;j) {xj - Xkf (1.5) 

j=l l<j<k<N 

where / = (a, 0), — oo < a < < +oo and a;(a;) is a weight function on the interval 
/ all finite moments of which exist. 

Remark 1.1 (1.5) was imposed on a definition of a matrix ensemble in [Me], and 
obtained in [TW] for 

F{A) = exp( - Tr{V{A))) (1.6) 

where V{x) is a real-valued function such that uj{x) = exp(— T^(a;)) defines a weight 
function. In addition, it is easy to draw the conclusion similar to (1.5) for the 
orthogonal and symplectic ensembles. 

Remark 1.2 The appearance of the characteristic function xi of ths set / means 
that we just consider sub-ensemble of Hermitian matrices space whose eigenvalues 
are all in /. Moreover, I can be given by the union of some disjoint intervals and 
then Lo{x) is the associated weight function. 

Prom (1.5) the n— point correlation function is defined in [Me] by 

Rn{xi, ...,Xn)= [ ■■■ [ Pn{xi,X2,...,Xn) dXn+l...dXN- (1.7) 

{N-n)\Jj J I 
If we introduce the orthogonal polynomials 

Pj{x)pk{x)uj{x)dx = 6jk, j, fc = 0, 1,... (1.8) 

and the associated functions 



fkix) =_pfe(x)v/t^(x), (1.9) 

then by the property of Vandermonde determinant the joint distribution function 
(1.5) reads: 

Pjv(a;i,a;2,...,a;jv) = -^(det[<^j_i(a;fe)]|^fc=i)^ 
1 

m' 

where 

KN{x,y) = ^ ^pj{x)^pj{y). (1.11) 

3=0 



= —Aei[KM{x,,Xk)]\lk^, (1.10) 
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Because of the orthonormality of the (pj{x)'s one can show just as in [Me, Chap. 6] 
that 

Rn{xi,...,Xn) =dct[KN{Xj,Xk)]\lk=l- (1-12) 

It is of considerable interest to obtain the behavior of n— point correlation of 
(1.12) after some appropriate scaling in the limit of large N . In particular, putting 
n = 1, we get the density of eigenvalues (also called level density ) 

N-l 

Ri{x) = Kn{x,x) = ^ ^^j{x) (1.13) 

and denote the normalized density by 

1 1 

nM^) = nE^M (1-14) 



Remark 1.3 For the eigenvalues xi,X2, - ■■ ,Xn of a random Hermitian matrix 
A, the integral of the normalized counting function over an interval A = (c, d) is 
given by 




(1.15) 



A calculation in [PS] shows that 

E{un){A) = ^Ri{x)dx (1.16) 

where E{un){A) denotes the expectation of z>'jv(A) with respect to the probability 
measure of (1.5). 

Our motivation is to obtain the density in the scaling limit of large N from 
(1.14) and further study its polynomial-perturbation invariance. Before our results 
are stated we review some known results about the density. 

As we know, Wigner in [Wigl, 2] not only got his famous semicircle law (here 
corresponding to the weight u}{x) = called Gauss unitary ensemble, denoting 
GUE) 

aix) = -Vl-x^ X[-i,i]ix), (1.17) 

but also invented the calculation method which had some independent interest as 
he thought. Afterwards, Bronk got the associated density for lu{x) = x^e~'^, p > 
— 1,0 < X < +00 for Laguerre ensembles in [Br] and Leff for u){x) = (1 — a;)"(l + 
x)^, a,b > —1, —1 < a; < 1 for Jacobian ensembles in [Lcf]. Recently, an irradiativc 
idea in [HT] was introduced by Haagerup and Thorbjornsen to give a short proof of 
Wigner's semicircle law, using the Laplace transform of (1.14) and the property of 
Hermitian polynomials. In [Lef], Lcdoux pushed forward the investigation in [HT] 
and in a slightly different way obtained Wigner's semicircle law, also the densities 
for Laguerre and Jacobian unitary ensembles but whose expressions appeared com- 
plicated. In addition, based on the spirit of statistical mechanics the method of 
equilibrium measure is used to obtain the density, mainly for the weight function 
w{x) = ex.p{—V{x)), 

Vix) = J2mx''"' + • • • + 70, 72m > 0. (1.18) 
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See [Jo] , [DMK] or [De] for the equilibrium measure method. In a recent survey on 
orthogonal polynomial ensembles ([K6]), Konig obtained semicircle law for GUE 
respectively by moment method and equilibrium measure method. 

In this paper, we will generally deal with the weight function u!{x) on the inter- 
val / and obtain the density of eigenvalues (t{x). In fact, Professor W. Van Assche 
told us that P. Nevai and W. Van Assche had shown that the density and zero 
distribution of the orthogonal polynomials were the same ([Val], Theorem 5.3), see 
Remark 1.8 below. However, our method is to calculate directly the moments with- 
out other knowledge than three-term recurrence formula , and further the density 
was obtained using a different method. To our knowledge, the family of densities 
satisfying differential equations of first order is first obtained in the present paper. 

It is a well known result that three- term recurrence formula holds for the orthog- 
onal polynomials defined by (1.8) 

xpnix) = an+ip„+i{x) + bnPnix) + a„p„_i(a;), n = 0, 1, • • • (1-19) 

where o„ > 0,p-i{x) = 0. We will assume that there is a positive and non- 
decreasing sequence c„ such that 

lim — = a > 0, lim — = b. (1.20) 

n— ►oo Cn n^oo c„ 

An extra condition on the contraction sequence is that c„ is a regularly varying 
sequence with index A > 0, i.e. 

c„ = n^L{n) (1.21) 
where L : (0, -|-oo) — > (0, -|-oo) is slowly varying , that is, 

lim = 1, \/t>0. (1.22) 

x^oo L(x) 

The condition of (1.21) was first introduced by W. Van Assche to study the asymp- 
totics for orthogonal polynomials, see [Va2] as a general survey for the condition of 
(1.20) and [Fe, Chap.VIII] for regular functions (1.22). 

Remark 1.4 Assuming c{x) is a positive, non-decreasing differentiable function 
on the interval (0, -|-oo) with 

lim ^^ = A>0, (1.23) 

x^+ca c[x) 

then one easily knows that c{x) can be represented as 

c{x) = x^L{x). (1.24) 
A discrete version of this result is that (1.21) holds if 

lim nf^- l) =A>0. (1.25) 



Remark 1.5 The conditions of (1.20) and (1.21) are connected close with the 

asymptotic problems of the orthogonal polynomials. The class of Freud weights 
plays a most close role on (1.21), for examples, lo{x^ = exp{—Q{x)) where Q{x) 
grows like a power at infinity, in particular, 

Q{x) = 72ma;2™ + • • • + 70 , 72m > (1.26) 
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with A = 1/(2to), see [DKMVZ] ; 

uj(x) = |x|^e-l^l°, /3 > -1, a > (1.27) 
with A = 1/a, see [LMS] or [Val]. In addition, A = 1 for Laguerre weights 

uj{x) = a;"e-"=, a > -1 (1.28) 

and A = for Jacobi weights 

uj{x) = (1 - .t)"(1 + xf, a,f]> -1. (1.29) 

A classic result of E.A. Rakhmanov in [Si] asserts that if the weight ui on [—1,1] 
satisfying a; > a.e. on [—1, 1], then u){x) belongs to the Nevai-Blumenthal class, 
that is 

lim an = 1/2. lim b,, = 0. (1.30) 

n — >oo n — ^OC' 

Obviously, A = for the Nevai-Blumenthal class (sec [Va2]). We strongly refer the 
reader to [Lub] for a recent survey for a wide variety of weights on finite or infinite 
intervals. 

Now we can state our main results. In Section 2 we will introduce ascending, 
equilibrating and descending operators which describe the transforming of polyno- 
mials, and explicitly calculate the moments of the density. Then in Section 3 we 
consider a simple deformation of the moments for any given density and obtain a 
corresponding density determined by a differential equation with respect to the new 
moments. In section 4 using the results in Section 3 we give the proofs of Theorems 
1.2 and 1.5 below. 

First let us rescale the density of (1.14) by 

<^n{x) = ^Ri{cnx). (1.31) 
Note that a^ix) is our main object and we will study its limit behavior. 



Theorem 1.1 Denote the fcth moment of the scaling density crjv(a;) by m'^\ 
Under the contraction conditions of (1.20) and (1.21), we have 

lim M^^^ =Mk, fc = 0,l,... (1.32) 

iV— >(x> 

where 

, [fe/2] 

= TTAfc( ^ CiCi_^a^^h^-^^). (1.33) 

Remark 1.6 One sets a = 1/2 in the following since c„ can be chosen freely from 
some constant. Observe 

^ qq_^.a2^-6'=-2^- =L,(az+°-+ h) (1.34) 
j=o 

where the operator Lo{f) represents the constant term of Laurent series, using 
Cauchy contour integral, that is 

= 2^ / ^'-''^ 
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Theorem 1.2 A probability density cr(x) with its kth moment Mk exists, and is 
uniquely determined by the following differential equation of first order 

aix) - X[xa{x)f^ = ^ ^ Xu, (1-36) 

y/l — {x — by 

with the following conditions 

/+00 
a{x)dx = 1 (1.37) 
-oo 

where lb = (— 1 + &, 1 + &) and xj,, a characteristic function of 7^. Exactly, the 
support of a{x) can be restricted to a finite interval, that is, for A = 

^W = ^ ] (1-38) 

TT y'l - (a; - 6)^ 

while for A > 

supp(c7) = [Bi,B2] (1.39) 

where 

Bi = mm{b - 1, 0}, B2 = max{6 + 1, 0}. (1.40) 
One directly solves the equation of (1.36) and easily obtains 

Corollary 1.3 For 6 = and m € N, we have 

(1) for A = l/2m 

and 

(2) for A = l/(2m- 1) 

^ ' 2tt V2' \x\ 2 iCL 2' J ^ ' 

j=l J 2j 



For 6 = -1 and g G N U {0}, we have 
(3) for A = l/{q + 1) 



and 

(4)for A = l/(g+i) 

Remark 1.7 For the weight ij{x) = exp(— ^(a;)) from (1.18) which is correspond- 
ing with A = 1 /2to, the density was given in [Jo] in the form of r{x) {x — xi){x2 — x?). 
Here xi < X2 and r{x) is a polynomial of degree 2m — 2 depending on V{x)] when 
A = 1 and 6 = 0, it was obtained in [CI] using Mathematica for Meixner-Pollaczek 
polynomials. It is just the case where A = 1 and h = —1 for Laguerre weights of 
(1.28), see [Sze]. The authors believe that the densities in Corollary 1.3, especially 
for small integrals m and g, may appear in some other random matrix models. 
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From Theorems 1.1 and 1.2, it is obvious that 

Theorem 1.4 Denote (Jn{x) and a{x) as above , then 

crjv(a;) a{x) (1-45) 
where W means in the weak sense. 

Remark 1.8 Let xi^n,X2^m ■ ■ ■ ,Xn,n be zeros of p„. P. G. Nevai and J. S. De- 

hcsa in [ND] also got the same moments of zero distribution under the contraction 
condition of (1-20) (this fact was pointed out by W. Van Assche in [Va2]), namely, 

lim i^(^)'= = Mfe. (1.46) 

n— ►oo n ; Cn 

A beautiful probabilistic interpretation in [Val] shows Mk is the fcth moment of 
Nevai-Ullman measure. In addition, Nevai and Van Assche proved that weak con- 
vergence of zero distribution implied some weak convergence to the same probability 

Af-l 

measure for Christoffel functions, which are defined by (^ p'^{x))^^ (see Theorem 

i=o 

5.3 in [Val]). It is obvious that the density of eigenvalues and zero distribution are 
the same according to (1-46) and Theorem 1.1, thus we give a new proof. 

Let p(x) be a fixed Zth order polynomial, and we consider a new weight function 

Luix) ^p^ix)oj{x) (1.47) 

with 

Pj{x)pk{x)cd{x)dx = Sjk, j, fc = 0, 1, (1-48) 



Associated 1-point correlation function is given by 

N-l 

R^{x) = J2f{x)p^{x)Lu{x). (1.49) 

j=0 

Under the same scaling we write 

^n{x) = ^Ri{cnx). (1.50) 
Calculate the moments of (1.50) for the new weight function, we find 

Theorem 1.5 Denote the fcth moments of (T7v(.x) and (7n{x) by m\^^ and M^^\ 
respectively. Under the contraction conditions of (1.20) and (1-21), we have 

lim Mi^'> = lim M^^^ = Mk, A; = 0, 1, ... . (1.51) 

n— >C!0 n— >C!0 

Namely, 

^n{x) a{x) (1.52) 
where a{x) is the density determined by (1.36) and (1.37). 
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2. Calculation of the moments 

By the recursion formula of (1.19), we regard the multiplication by x as an 
operator A^, and it can be represented as 

A^=A+ + Ao + A_ (2.1) 

where A+,Ao and A^ are called ascending, equilibrating and descending operators 

respectively, defined by 

A+pn{x) = a„+ipn+i{x), AoPn{x) = bnPnix), ^_p„ (x) = a„p„_ i (x) . (2.2) 
Thus we calculate the A:th moment of aN{x) using (1.31), (2.1) and (2.2) as follows: 



M^^' = j x''(7N{x)dx 



1 f 

^ N-l 

' j:{i^^ + ^o + A.)^p„p,). (2.3) 



Let be a set composed of those terms in the expansion of (^+ + + AS)^ , 
in which the operators and A- exactly appear q times respectively. Note that 
(Tpj,Pj) j^2(^^-^ = for T ^ (J A^, then one obtains the following lemma: 



Lemma 2.1 

. N-l [fe/2] 

For convenience of the following calculation, we introduce one property of the 
regular varying functions. 

Lemma 2.2 Let c(x) be a positive and non-decreasing function on (0, +00), and 

c{x) = x^L{x) (2.5) 



where L{x) satisfies 



Then 



lim ^ = 1, Vt>0. (2.6) 

;— >+oo L(x) 



n^oo (c(n))*= 1 + Afc ^ ' 



Proof By Lebesgue's dominated theorem and exchanging limits and integrals , it 
is easy to prove using (2.5) and (2.6). □ 
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Now we give a proof of Theorem 1.1 by two cases of b. 
Case 1 b^O. 

Assume that 6 > for convenience. Write 

^ = a(l + e„), ^ = b{l + vn). (2.8) 

If one writes for j > k 

«i = . max {\im\, \Vm\}, (2.9) 

j — k<'m<o+k 

then 

lim Uj = 0. (2.10) 

j-»oo 

Thus we can suppose Uj < IjVj = 0, 1, . . . and by the definition of Uj one obtains 
max {am} <a Cj+k{^+Uj), max {bm} <b Cj+k{^ + Uj). (2.11) 

j—k<m<j+k j—k<m<j+k 

Furthermore, for T e A| we have 

(Tpj,pj) < ( max {a,n}) ( max {bm,}] 

\ / L^{w) \j-k<m<j+k J \j-k<m<j+k ) 

< a^%''-^''icj+k)\l+u,)\ (2.12) 
Summing by q and j, we get 

JV-l 

[k/2] jr E (cj+fc)'=(i + wj)'' 

Mr < ( E GlCl_y%>'-^'^) . (2.13) 



By Cauchy-Maclaurin summation formula and Lemma 2.2 one obtains 



E'(c,+fc)'=(i + Uj)'' i E\cj+k)' 



lim -7 = Hm , 

JV^oo [Cn) N^oc \Cm) 

N 

lim 



N-,oo {c{N)Y 
1 



1 + Afc 



(2.14) 



Analogously, for j > k, one obtains 

[fe/2] 

> ( E ^.'C'L,a^^6'=-^) (2.15) 



,[fe/2] ^ E (c,-fc)''(l-W,)* 
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and 



lim 

N^oo 



N-1 
j=k 



1 



Combining (2.13) — (2.16) we complete the proof of this case. □ 
Case 2 6 = 0. 
Similarly, write 

— = a(l + ?n), — = r]n 



(2.16) 



and 
then 



Vj ^ . max \rjm\}, 

j — k<m,<j-\-k 



lim =0. 



(2.17) 
(2.18) 
(2.19) 



Thus we can suppose < 1, Vj = 0, 1, For Q<j — k<m<j + k, one obtains 

o Cj-k{l - Vj) <am<a Cj+fe(l + vj), \bm\ < Cj+kVj. (2.20) 



Now we give an estimation of Mjf'^ when k is odd. For T G A?, by (2.20) we 



have 



If denoting 



then 



[fe/2] 



k~2q 



9=0 



lim Vj = 0. 



Thus by the Stoltz formula one obtains 



(iV) 



< 



JV-l 
j=0 



0. 



When k is even. For T G A^, g 7^ k/2, 



and for T e A^^^j > k, 



a'(ci-fc)'(l - < {TPhPj),,, , < aHcj+k)H^+Vjf 



If one writes 



Mi") 



N-l 



E E E (^^'. 



(2.21) 

(2.22) 
(2.23) 

(2.24) 

(2.25) 

(2.26) 
(2.27) 
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and 

N-1 



then 



AfW=Mf^+ii?]^^\ (2.29) 



Similar to the case where k is odd, by (2.25) one obtains 



(2.30) 



and similar to the case where 6 ^ 0, by (2.26) one obtains 

= T^^k"^'- (2-31) 

Combining (2.24), (2.29), (2.30) and (2.31) we complete the proof of this case 
where 6 = 0. □ 

Therefore, Theorem 1.1 has been proved. 

3. Deformation of the moments 

Following Wigner's original papers [Wigl] and [Wig2], to obtain the distribution 
of eigenvalues in random matrix models, a standard procedure is first to calculate 
the moments of all orders after some appropriate scaling of eigenvalues and then 
to determine an explicit distribution function of the moments. In this section we 
discuss some deformations of the moments of a given density function and determine 
the corresponding density function with respect to the new moments, and then 
give an application in the following section. A natural question: which kind of 
deformations of the moments has a corresponding distribution function? We will 
make an interesting try. 

Let f{x) be a continuous density function on the interval / = (a, /3), — oo < 
a < (3 < +00, and its moments 

nik = J^x''f{x)dx, fc = 0,1,2,... (3.1) 

exist. Now consider the following deformation of moments with a parameter A > 0, 

Mk = rT-rr'^k, k = 0,1,2,.... (3.2) 

Question: When is there a unique density function cr{x) which is a solution of 
the moment problem of (3.2), i.e. 

f+OO 

x''a{x)dx = Mk, k = 0,1,2,..., (3.3) 



/ 



and further how to determine cr(x) from the given density f{x)l 

Note that two density functions f{x) ~ g(x) are said to be equal if their distri- 

f{s)ds and / g{s)ds are equal at their points of continuity. 

-oo J — oo 

To make sure that the density is unique we assume the moments ruk, k = 0, 1, 2, . . . 
satisfy Carleman's condition 



E-2r' = oo. (3.4) 

k=0 
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It is obvious that the moments Mk also satisfy Carleman's condition 

oo 

^M-V^'= = oo. (3.5) 

k=0 

Due to Carleman's famous theorem (see [ST] or [Fe]) (3.5) assures that the moment 
problem of (3.2) is determined, i.e. numbers Mk determine a unique density whose 
kth moment is M^- 

To prove that the moment problem of (3.2) has a solution a{x) whose spectrum 

support supper (x) is to be contained in the interval J, given in advance, we first 
introduce an important theorem. Let P{u) be any polynomial in u, 

P(u) = ^xfe«'= (3.6) 
fe 

where numbers Xk are real constants. Introduce the functional /u(P) defined by 

tx{P) = Y,MkXk- (3.7) 

k 



Theorem 3.1 ([ST], Theorem 1.1) A necessary and sufficient condition that the 
moment problem defined by the sequence of moments Mfc shall have a solution on 
J is that the functional ^(P) be non-negative, that is 

At(-P) > 0, whenever P{u) > on J. (3.8) 

Note that Theorem 3.1 can be applied to derive explicit necessary and sufficient 

conditions by a special choice of J, which depend on representations of non-negative 
polynomials on J. In particular, it is a well-known fact (see [P6S]) any polynomial 
P{u) > for all real u can be presented as 

P{u) = P,{uf + P2iuf (3.9) 

where Pi (w) and Pj (u) are polynomials with real coefficients. 

If we take for P{u) the particular polynomial P{u) = {xq + xiu -\ h XnU")"^, 

we have 

n 

fi{P) = J2 Mj+kXjXk. (3.10) 

j,k=0 

From the theory of quadratic forms and (3.9) it is well known that the conditions 
of (3.8) are equivalent to 

\M,+k\lk=o>^, n = 0,l,2,-.. (3.11) 
if the spectrum of the solution is not reducible to a finite set of points. 

Now we can state our theorem as follows: 



Theorem 3.2 A probability density a{x?) on the real axis with its fcth moment 

Mk of (3.2) exists, and is uniquely determined by the differential equation of first 
order 

a(,T)-A[xa(x)]« =/(x)x/ (3.12) 
with the following conditions 

a{x)dx = 1. (3.13) 

-OO 
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Proof. It is sufBcient to prove (3.11) if {Mj+k)J^k=o ^ positive definite matrix. 
Taking f{x)xi for the solution of the moments ruk on the real axis, we have 

\mj+k\lk=o > 0, n = 0,1,2,- •• , (3.14) 

that is, 

A„ = {mj+k)^k=o is a positive definite matrix, n = 0, 1, 2, • • • . (3.15) 
The formula of (3.2) shows that 

{Mj+k)l„^0=An*An (3.16) 
where * represents Schur product and 

^" = (tTa^)".-o> n = 0,l,2,.... (3.17) 

Thus by the property of Schur product it is sufficient to prove that A„ is a positive 
definite matrix. Note that 

ITS = 

then one can obtain 

"^^^^^ = T^TTw /'•••/' n - n > 0. (3.19) 

To derive the equation of (3.12) we make a Fourier transform and write 

H{t) = r e'*''a{x)dx, F{t) = f e'^''f{x)dx. (3.20) 

J —OO J I 



Thus, 



/ itx e^^='a{x)dx = ^ ^ ' Mk+i 



k=0 

fe+1 



k=0 ^ ' 



E^^^'^- (3-21) 



Obviously, 



/ + 00 OO /-.U: 

(1 + Aite) e'*V(x)dx = ^^(l + AA;)Mfc 
-°° fc=0 



fe=0 



^ e'*''f{x)dx. (3.22) 



Namely, 

H{t) + XtH '{t) = F{t). (3.23) 
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We also make an inverse Fourier transform. Combining 

r"+°° 1 r -, (1) 



— e-^^''tH' {t)dt = -— / e-''''H '{t)dt 

/+00 
e-it^H{t)dt 
-OO 



(1) 



1 

= -[xa{x)]^^'' (3.24) 
and (3.23) one obtains the differential equation on real axis 

aix)-X[xaix)]^^'^ = f{x)xi (3.25) 
whose possible singular points are 0, a and /? (if a and (3 are finite numbers). □ 

Note that if a or /? is finite we have an exact information about the spectrum 
support J of a{x). 

Case 1. — OO < a < 0, (3 = +oo (or a = — oo, < /3 < +oo similarly ). J = I. 

Any non-negative polynomial P{u) on / = (a, +oo) can be represented by 

P{u) = Pi{uf + P2{uf + [u - a){P3{uf + Pi{uf) (3.26) 

where Pi(w), -P2('«), Pd.{u) and Pi{u) are polynomials with real coefficients. 

An analogous procedure to the above arguments of (3.9) — (3.11) one obtains 
the conditions of (3.8) are equivalent to the matrices 

{Mj+k)lk=o^ n = 0,l,2,..- (3.27) 

and 

{Mj+k+i-aMj+k)lk=o. n = 0,l,2,--- (3.28) 

are positive definite . 
Note that 

- "^'^ = 1 + A(fc + 1) + ^-^" \l + Afc)(l + A(fc + l)) 

and 

{m,+k)lk=o^ {mj+k+i - amj+k)lk=o^ n = 0, 1, 2, • • • (3.30) 
are positive definite, by the property of Schur product one obtains the matrices of 
(3.27) and (3.28) are positive definite. 

Case 2. -oo <a<0,0</3< +oo (or -oo < a < 0,Q < /3 < +oo similarly ). 
J = I. 

Any non- negative polynomial P{u) on I = {a, (3) can be represented by 

P{u) = Pi{uf + {u- a){(3 - u)P2{uf (3.31) 

where Pi{u) and P2{u) arc polynomials with real coefficients. 

In this case the conditions of (3.8) arc equivalent to the matrices 

iM,+kyik=o, n = 0,1,2,- •• (3.32) 

and 

{-Mj+k+2 + (/? - a)Mj+k+i - apMj+k)lk=o, n = 0, 1, 2, • • • (3.33) 
are positive definite . 
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Note that —Mk+2 + {P — a)Mk+i — aj3Mk can be rewritten by 
-mfe+2 + (/? - a)mk+i - a/3mfe _^ [-a(3) 



1 + A(fc + 1) ' ' ^'(1 + Afc)(l + A(fc + 1)) 

(l + A(A: + l))(l + A(fc + 2)) ^ • ' 

and 

{mj+k)'j,k=o^{-'n^k+2 + {l3-a)mk+i-a(3mk)'-^k=a^ n = 0, 1,2,--- (3.35) 

are positive definite, again by the property of Schur product one obtains the ma- 
trices of (3.32) and (3.33) are positive definite. 

Remark 3.1 If ajS > 0, we cannot make sure that the matrices of (3.33) are 
positive, e.g., a = 1, /S = 2,n — 

-M2 + (/3 - a)Mi - aPMo = + - 2mo < (3.36) 

for the sufficiently large A. 

Case 3. 0<a</J(ora</3<0 similarly ). J = (0, /?) (or J = (a, 0)). 
We take f{x) for the density function on (0,/3) 



f{x) a <x < P 

< a; < a (3.37) 



and by using Case 1 and Case 2 it is obvious. 

Anyway, we get 

Theorem 3.3 A probability density (j{x) with its fcth moment Mk of (3.2) exists, 
and is uniquely determined by the following differential equation 

(1) a < < /3, 

a{x) - A[a;c7(a;)](^) = f{x), x G I; (3.38) 

(2) a < /3 < 0, 

a{x)-X[xaix)]'^'^ = f{x)xi, x G {a,Oy, (3.39) 

(3) < a < /?, 

a{x) - A[a;a(x)](i) = f{x)xi, x € (0,/3). (3.40) 

4. Proofs of Theorems 1.2 and 1.5 
Proof of Theorem 1.2. 

In [Wigl] Wigner made use of the integral representation of Bessel function of 

order 1 (pointed out by W.FcUcr to him) to get his semicircle law. However, in the 
original analysis he got the semicircle law by leading a differential equation. In the 
following we will combine these two kinds of method to derive the density. 
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k 

Uniqueness: Note that using (1.33) dLcik-i < ^'^'^ Isolds. Writing B = 

j=o 

a+\b\, then we have 

k 

< E ciA-j^'^b''-'' < m"'- (4-1) 

Thus, the Carleman's condition is satisfied and the density function is determined 

by the moments. 

Derivation of the differential equation: Putting 2a = 1, using the integral repre- 
sentation of Bessel function of order zero (see [Sze] ) and calculating directly 



/+00 /-.U 
(l + Xitx) e**"(T(a;)dx = ^ i-j-(l + Afc)M, 



k=0 



fc=0 j=0 

\k 



E(itf ^ ( a \k 



fe=o 



fe=0 



= Z/0 exp {it{az + - + 6)^ 
= e'*^Lo exp {it{az+-^ 



__ itb \ " i'i-tay 

Mb 



= e'" Jo{2at) 



1 /■! pita; 

= e'*^i / -A dx 



= / e''''ft,ix)dx (4.2) 

J6-1 



where 



- (x - 6)2 

By using Theorems 3.2 and 3.3 in Section 3 one obtains 

a{x) - X[xc7{x)f^ = Mx)xi, (4.4) 
and supp(cr) = [Bi,B2] for A > while for A = 0, obviously 
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Now wc give an exact solution of the equation of (4.4). Note that when A > 
the equation of (4.4) is known as a Cauchy-Euler equation of order 1 (see [GN], 
P99). Thus we have 

(1) — 1 < 6 < 1. is a singular point of the equation of (4.4), and we have to 
determine how the solutions for a; < and a; > can be pieced together to give 
solutions valid on the whole interval /(, (see [GN], P22 ). So one obtains 

{1 1 f''^^ 1 
-.T^"^ / s^i fb{s)ds x>0, 
-{-Xp-^J^ 1 (-S)-V6(s)rfs X<0. 

rb+1 

Note that xa{x) is absolutely continuous on lb. Thus / a{x)dx = 1 from the 

Jb-l 

equation of (4.4) by using a{b — 1) = <j{b + 1) = 0. Besides, a{x) is continuous on 
lb when < A < 1, while is a singular point of a{x) when A > 1. 

(2) 5 = ±1. 
If & = 1, then 

a{x) = \x^-^ j s-i J. ds, a;e(0,2). (4.7) 

While 6 = —1, we have 

a(x) = \{-x)\-^ r {-s)-i ] ds, xe(-2,0). (4.8) 

A J -2 ^y-s{s + 2) 

(3) 6 > 1. The equation of (4.4) can be rewritten by one nonhomogeneous linear 
equation on /ft = (6 — 1, b+ 1) 

a(x) - X\xa(x)] = ^ (4.9) 

and the other homogeneous linear equation on (0, 6—1) 

a{x) - X[xa{x)]^^'^ = 0. (4.10) 
Next we solve the related nonhomogeneous linear equation of (4.9) 

a{x) = -x^-^ j s-ifbis)ds (4.11) 

and the related homogeneous linear equation of (4.10) 

a{x) = C+x^-^ (4.12) 

where 

C+ = - s-ifb{s)ds. (4.13) 

^ Jb-i 

Note that C+ assures that a{x) is continuous on (0, 6+1) and J^^^ a{x)dx = 1. 

(4) 6 < — 1. Similar to the case where 6 > 1. One obtains the solutions of the 
equation of (4.4) on /b = (6 — 1, 6+1) 

a{x) = \{-x)-^-' r {-s)--^fb{s)ds (4.14) 
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and on (6 + 1, 0) 

(j{x) = C-{-x)i-'^ (4.15) 

where 



1 



b+i 



C- = - (s)-^ f,{s)ds. (4.16) 

^ Jb-l 

In the end we will complete the proof of Thcorcni 1.5. First, setting 

Hn = span{pQix),pi{x), ■ ■ ■ ,p„_i(a;)}, (4.17) 

then Hn is an n-dimensional subspace of L'^{uj). It is obvious that po{x)p{x), • • • , 
Pn-i-i{x)p{x) is a family of normalized orthogonal vectors in Tin, extended by 

4"'(^)>--- ,e'il\{x),po{x)p{x),--- ,p„-i-i{x)p{x) (4.18) 

to a normalized orthogonal base of Hn ■ 

Let P„ be a projective operator from i^(w) to Hn- We construct an operator 
from Hn to itself as follows, 

T^'^) =PnoA^^:Hn^Hn (4.19) 
where is the multiplication by x. 

Lemma 4.1 Denote the A;th moments of aN{x) and aN{x) by Mjf^^ and M^^^ 
respectively, then 

Mf'=Mf) + e^ (4.20) 

where 

^ JV-l l-l 

®-=Ar7b^( <^x(ft^>),ftP)..(.)-E(^l^'^(4"^4"\^(.))- (4-21) 

j=N-l j=0 

Proof. We first point out that 

M^^' = J x''aN{x)dx 



On the other hand, by the normalized orthogonal base of (4.18) one obtains 



(4.22) 



N-l-l l-l 

= E {T'N\hp\hp) + E {T'N\ef\ ef )),,(^,. (4.23) 
j=o j=o 
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Thus, 



— ^fe ^ I x''fAx)p^{x)uj{x)dx 



Af-/-1 ^ l-l 



^( E {j^x''fA^)p\x)uj{x)dx + J^J^x\e^P{x)fuj{x)dx)+eN 
1 

N 

= M;^"^^ + e;v. (4.24) 



N (cn) - 

JV-i-l l-l 
j=0 j=0 



□ 



Lemma 4.2 Write || • || = (•, and for = 1, 2, • • • , 

k-l 



where 



P^/ll<3'=(n^"+0 11/11' (4.25) 



max{a,-, (4.26) 

0<j<n 



N-1 

Proof. Set f{x)= Yl hPji^) and note that 

JV-l 



J=0 
AT-l 



= II E + + '^J^'i-i 

= II E*^^J-1"J + ^i^i + 'j+i«j+ibill 

= E i^j-i^i + ^j^^' + ^j+i^j+ii^ 



< ^DlWff- (4.27) 
Thus, IjA^/ll < 3D„||/|| and (4.25) is easily proved. □ 

Proof of Theorem 1.5: By (1.21) and Lemma 4.2, there exist constants C\ and 
C2 which only depend on A; such that 

|(^S(P,P),P,-P)^.(„)| < \\Al{pjp)\\ < 3^i)^+,+,_2 < C,{cNf (4.28) 
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and 

|(rf (ef )),ef ))^,(^)| < IIP^II llA^ef )|| < 3'=I?^+,_i < C^icM)'. (4.29) 
Thus, using (4.21), for the large A'', 

^ ' j=N-l j=0 

Again by Lemma 4.1 we get 

lim Mi^'> = Hm m'^' = M^. (4.31) 

□ 
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